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Abstract 

We calculate the amplitude for the forward electroproduction of two light vector mesons 
in next-to-leading order BFKL. This amplitude is written as a convolution of two impact 
factors for the virtual photon to light vector meson transition with the BFKL Green's 
function. It represents the first next-to-leading order amplitude ever calculated for a 
collision process between strongly interacting colorless particles. 



1 Introduction 



Collision processes between strongly interacting particles in the limit of large center-of- 
mass energy are the best candidates for the realization of the BFKL dynamics [1], provided 
that a "hard" scale exist which allows the use of perturbation theory in the strong coupling 
ots- 

In the BFKL approach, both in the leading logarithmic approximation (LLA), which 
means resummation of leading energy logarithms, all terms (Q;sln(s))", and in the next- 
to-leading approximation (NLA), which means resummation of all terms ^^(q;^ ln(s))", 
the (imaginary part of the) amplitude for a large-s hard collision process can be written 
as the convolution of the Green's function of two interacting Reggeized gluons with the 
impact factors of the colliding particles (see, for example. Fig. 1). 

The Green's function is determined through the BFKL equation. The NLA singlet 
kernel of the BFKL equation has been achieved in the forward case [2], after the long 
program of calculation of the NLA corrections [3] (for a review, see Ref. [4]). For the non- 
forward case the ingredients to the NLA BFKL kernel are known since a few years for the 
color octet representation in the t-channel [5] . This color representation is very important 
for the check of consistency of the s-channel unitarity with the gluon Reggeization, i.e. 
for the "bootstrap" [6]. Recently the last missing piece has been determined for the 
determination of the non-forward NLA BFKL kernel in the singlet color representation, 
i.e. in the Pomeron channel, relevant for physical applications [7]. 

On the other side, NLA impact factors have been calculated for colliding partons [8] 
and for forward jet production [9]. Among the impact factors for transitions between 
colorless objects, the most important one from the phenomenological point of view is 
certainly the impact factor for the virtual photon to virtual photon transition, i.e. the 
7* 7* impact factor, since it would open the way to predictions of the 7*7* total 
cross section. Its calculation is rather complicated and only after year-long efforts it is 
approaching completion [10]. 

A considerable simplification can be gained if one considers instead the impact factor 
for the transition from a virtual photon 7* to a light neutral vector meson V = p^,u}, (p. 
In this case, indeed, a close analytical expression can be achieved in the NLA, up to 
contributions suppressed as inverse powers of the photon virtuahty [11]. In particular, 
it turns out that (a) the dominant hclicity amplitude is that for the transition from 
longitudinally polarized virtual photon to longitudinally polarized vector meson; (b) the 
impact factor, both in the LLA and in the NLA, factorizes into the convolution of a hard 
scattering amplitude, calculable in perturbative QCD, and a meson twist-2 distribution 
amplitude [11]. 

The knowledge of the 7* — > impact factor allows for the first time to determine 
completely within perturbative QCD and with NLA accuracy the amplitude of a physical 
process, the 7*7* — > VV reaction. This possibihty is interesting first of all for theoretical 
reasons, since it could shed light on the role and the optimal choice of the energy scales 
entering the BFKL approach. Moreover, it could be used as a test-ground for comparisons 
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with approaches different from BFKL, such as DGLAP, and with possible next-to-leading 
order extensions of phenomenological models, such as color dipole and /ct-factorization. 
But it could be interesting also for the possible applications to the phenomenology. Indeed, 
the calculation of the 7* — K impact factor is the first step towards the application of 
BFKL approach to the description of processes such as the vector meson electroproduction 
7*p — > Vp, being carried out at the HERA collider, and the production of two mesons in 
the photon collision, 7*7* — > VV or 7*7 — > VJ/^, which can be studied at high-energy 
e'^e" and 67 colliders. 

In this paper we concentrate on the NLA forward amplitude for the 7*7* VV 
reaction. Such a process has been studied recently in Ref. [12] in the Born (2-gluon 
exchange) limit for arbitrary transverse momentum. First of all, we show how the available 
results for the 7* — > ^ impact factor and the BFKL Green's function can be put together 
to build up the NLA amplitude of the 7*7* — > VV process in the MS scheme. Then we 
restrict ourselves to the particular case of collision of virtual photons with equal virtualities 
and present some numerical estimates of our result, aimed at showing the extent of the 
contributions to the NLA amplitude from the impact factor and from the NLA kernel 
and the dependence on the energy scale introduced in the BFKL approach and on the 
renormalization scale which appears in the MS scheme. 

Moreover, we show that, despite being the NLA corrections large and of opposite sign 
with respect to the leading order, it is possible to achieve a well-behaved form of the 
amplitude, by a suitable choice of the energy and renormalization scale parameters. The 
procedure we adopted here to optimize the perturbative result will be compared with 
other approaches to optimization and also with methods based on the improvement of 
the NLA BFKL Green's function [13] in forthcoming publications. 

When the present paper was ultimated, a new paper appeared [14] in which forward 
amphtude of 7*7* — > VV process is studied in the LLA and also an estimate for the NLA 
effects based on the specific assumptions is given. We found that at LLA level our results 
coincide with that ones in [14]. In order to clarify the physics behind the BFKL NLA 
corrections it is important to compare our exact NLA results with different approximate 
approaches to the NLA effects. However, such a study goes beyond the scope of the 
present paper. It will be given in a separate publication where, in particular, we plan to 
make a systematic comparison of exact NLA results with the estimates of [14]. 



2 The NLA amplitude 



We consider the production of two light vector mesons {V — p'^,uj,(f)) in the collision of 
two virtual photons, 

l*{p)l*{p')^V{pr)V{p^). (1) 
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Figure 1: Schematic representation of the amphtude for the 7*(p)7*(p') — > Vipi) V{P2) 
forward scattering. 

Here, pi and p2 are taken as Sudakov vectors satisfying Pi = P2 = ^ and 2{pip2) = s; the 
virtual photon momenta are instead 



Qi I I Q2 /o\ 
p = api p2 , p = ap2 —pi , (2) 



so that the photon virtuahties turn to be = —Ql and {p'Y = —Q\. We consider the 
kinematics when 

s > 2 > A|^^ , (3) 

and 

a = 1 + ^ + , a' = 1 + ^ + . (4) 

s s 

In this case vector mesons are produced by longitudinally polarized photons in the lon- 
gitudinally polarized state [11]. Other helicity amplitudes are power suppressed, with a 
suppression factor ~ mv/Qi,2- We will discuss here the amplitude of the forward scat- 
tering, i.e. when the transverse momenta of produced V mesons are zero or when the 
variable t = {pi — pY takes its maximal value to = —QiQl/s + 

The forward amplitude in the BFKL approach may be presented as follows 

J2 ^2 S+ico 

Im, (^) = / ^0) / 4|$2(-g2, .0) / ^ (-) GUqi, 92) . (5) 

(zvr)^ J qi J q2 J 2m \soJ 

5—ioo 

This representation for the amplitude is valid with NLA accuracy. Here $i(gi,so) and 
*^'2(— ^2, So) are the impact factors describing the transitions 7*(p) — > V{pi) and 7*(p') 
V{p2), respectively. The Green's function in (5) obeys the BFKL equation 

S'^iQi -q2) G^{qi, 52) - j d^q K{qi, q) G^{q, ^2) , (6) 
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where K{qi, ^2) is the BFKL kernel. The scale sq is artificial. It is introduced in the BFKL 
approach at the time to perform the Mellin transform from the s-space to the complex 
angular momentum plane and must disappear in the full expression for the amplitude at 
each fixed order of approximation. Using the result for the meson NLA impact factor 
such cancellation was demonstrated explicitly in Ref. [11] for the process in question. 

It is convenient to work in the transverse momentum representation, where "trans- 
verse" is related to the plane orthogonal to the vector mesons momenta. In this represen- 
tation, defined by 

QWi) = QiWi) , (7) 
(gllg-) = _ ) , {A\B) = {A\k){k\B) = I (fkA{k)B{k) , (8) 

the kernel of the operator K is 

K{q^,q^) = {q^\k\qr) (9) 
and the equation for the Green's function reads 

i = (u; - k)G^ , (10) 

its solution being 

G^^{u-k)-K (11) 

The kernel is given as an expansion in the strong coupling, 

k = ask^ + alk^ , (12) 

where 

OlgNc nn\ 

as = (13) 

TT 

and Nc is the number of colors. In Eq. (12) k^ is the BFKL kernel in the LLA, k^ 
represents the NLA correction. 

The impact factors are also presented as an expansion in ag 

where fv is the meson dimensional coupling constant {fp ^ 200 MeV) and Cq should be 
replaced by e/y/2, e/(3-\/2) and — e/3 for the case of p°, u and meson production, 
respectively. 

In the coUinear factorization approach the meson transition impact factor is given as 
a convolution of the hard scattering amplitude for the production of a coUinear quark- 
antiquark pair with the meson distribution amplitude (DA). The integration variable in 
this convolution is the fraction z of the meson momentum carried by the quark {z = 1 — z 
is the momentum fraction carried by the antiquark): 
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The NLA correction to the hard scattering amphtude, for a photon with virtuahty 
equal to Q^, is defined as follows 

1 /• , g2 



with t{z) given in the Eq. (75) of Ref. [11]. Ci^IIq'^) are given by the previous expression 
with replaced everywhere in the integrand by Qf and Ql, respectively. 

The distribution amplitude may be presented as an expansion in Gegenbauer polyno- 
mials 

(j)ll{z, i^f) = Qz{l - z) [l + a2{i^F)C2^\2z - 1) + a4i^F)Cl^\2z - 1) + . . .] . (17) 
The scale dependence of aninF) is well known [15]: 

aM = L^-'^'aM , (18) 

where L = as{iJiF) / ots^lA 

/5„^i^-f (19) 

is the leading coefficient of the QCD /9-function, with n/ the number of active quark 
fiavors. The anomalous dimensions 7„ are positive and grow with n. Therefore any DA 
approaches the asymptotic form (l)f{z) = 6z{l — z) at large /if-^ 

Below we will use the DA in the asymptotic form. Besides the simplicity of the 
following presentation, the reason is twofold. Presumably, the form of DA chosen at low 
Hf will affect mainly only the overall normalization of the amplitude but not the sum of 
BFKL energy logarithms and the resulting dependence of the amplitude on the energy in 
which we are primarily interested in this study. Another point is that, according to QCD 
sum rules estimates [16], 02(1 GeV) is 0.18 ± 0.10 for p and ± 0.1 for (f). Therefore 
may be indeed a good approximation for the DA of light vector mesons. Integrating over 
z in (15) with (j)\\{z, /ip) = (p^i^), we obtain, for photon virtuality Q^, 



(a = l^] =6a[l-^ln^ 



(20) 



where c = ^Ja-\-l/A . cf^ are given by the previous expression with Q"^ replaced by Q\ 

and respectively. For the NLA term Ci'^{q'^) the integration over z can be performed 
by a numerical calculation. 

To determine the amplitude with NLA accuracy we need an approximate solution of 
Eq. (11). With the required accuracy this solution is 



(21) 



"'^Thc dependence of the resulting amplitude on jip is subleading. Due to the collincar counterterm, see 
Eq. (72) of [11], the NLA correction to the meson impact factor contains a term proportional to ln(/iF), 
see Eq. (75) of [11], which compensates in the amplitude with NLA accuracy the effect of the meson DA 
variation with ^.p. 
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The basis of eigenfunctions of the LLA kernel, 



is given by the following set of functions: 



for which the orthonormality condition takes the form 



(22) 



(23) 



2^ 



6{u - v') 



(24) 



The action of the full NLA BFKL kernel on these functions may be expressed as follows: 



K\v) - «.(/^ii)x(i^)|i^)+«.'(/^i^)(x('H^) + :^xWln(/^^))k: 



dv , 



(25) 



where the first term represents the action of LLA kernel, while the second and the third 
ones stand for the diagonal and the non-diagonal parts of the NLA kernel. The function 
X^^\^), calculated in [2], is conveniently represented in the form 



f3o 



8 Nr. 



(26) 



where 



X(i^) 



^^X(^) - 6C(3) - x"(^) - 



+ 



2 / dx 



TT^ sinh(7rz/) 
2i/ cosh^(7ri/) 



cosh(7ri/ 
nf \ 11 + 12z/2 ^ 



m) 16(1 + ^^2) 



(i + x)V^ 



- Liofxl 



Liofx) 



X 

-I 



+ 4 0(1/) 



Here and below x'(^) = d{xi^)) / ^^'i x"{^) = d^ix{^))/d'^^- 



(27) 
(28) 



We will need also the \u) representation for the impact factors, which is defined by 
the following expressions 



^^^^ / d,.'c,{u')W\q), 

(J J 



+ 00 



duc2{v){^v), (29) 



cM = I d^qCr\q^)^^^ , c,{u) = I d^qC^^\q^)^^^^-j^ , 



(30) 



and by similar equations for c^\i') and C2\i^) from the NLA corrections to the impact 
factors, C[^\q^) and C^^\q'^). 

Using (21) and (25) one can derive, after some algebra, the following representation 
for the amplitude 



Irus (A) 
D1D2 



dz/(— J Oi^{i^R)ci{iy)c2{iy) 



1 + "^(/iij) 



+ 



, , 10 .dln{^) ■ 
-X(i^) + y + ^ + 2 Hl^l) 



We find that 



Ci{v)c2{v) = 



1 



^ r[3 ± 2iiy] cosh(7rz/) ' 



Q1Q2 \Ql J 32u{l + i/2) cosh3(7ri/) ' 



dln(^^ 
dv 



(31) 

(32) 
(33) 

(34) 



It can be useful to separate from the NLA correction to the impact factor the terms 
containing the dependence on sq and on /3o, 



CW(g2) 



dz — 



-»2 

q 



+ zzQ^ 



(35) 



X 



1 1 / -So \ , (ia + zz 



/So 



a'^z'^z^ 



Accordingly, one can write 

cSM = 5SM+cg(^), (36) 

where £1^2 (^) ^-r® the contributions from the terms isolated in the previous equation and 
c^li^) represent the rest. After straightforward calculations we found that 



\QiQ2r^ ' 2N, 



\Q1Q2 3 



+ -0(3 + 2iv) + -0(3 - 2iv) -il^l^ + iu^ - ^tjj - lu 



■ (37) 



Using Eq. (31) we construct the following representation for the amplitude 
Q1Q2 TrUsA 1 



D1D2 s 



(27r)^ 



X 



(38) 
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where the coefficients 



Q1Q2 



-|-oo 

j duci{u)c2{u) 



are determined by the kernel and the impact factors in LLA. Note that 

^'o = ^(7C(3)-6), 



(39) 



(40) 



therefore in the Born (the 2-gluon exchange) hmit our result coincides with that of 
Ref. [12]. 



The coefficients^ 



dn = n In 



So 



Q1Q2J ' 4iV, 



+ 



(n + l)^ln 
a 



bn-i , f I^R \ n{n- 1) 



Q1Q2 



dp (n+ l)f{v)ci{v)c2{u) 



Q1Q2J 2 



(n-iy. 



(41) 



Q1Q2 



diy Ci{u)c2{iy) 



X""^(^) 
(n- 1)! 



cj^V) -c^\u) x(^) 
c,{u) + C2{u) '>xiu) 



are determined by the NLA corrections to the kernel and to the impact factors. Here we 
use the notation 



f{jy) ^ + -0(3 + 2iiy) + -0(3 - 2iiy) - '4' [ - + - if^ 



(42) 



One should stress that both representations of the amplitude (38) and (31) are equiv- 
alent with NLA accuracy, since they differ only by next-to-NLA (NNLA) terms. Actually 
there exist infinitely many possibilities to write a NLA amplitude. For instance, another 
possibihty could be to exponentiate the bulk of the kernel NLA corrections 

XmsiA) s T /s\"=(''«M'')+"'('^«)(xM+JfexM[-xM+f]) , 



X 




1 + ^.(,.) ( ^ + ^ ) + ^li,,) m (A) ^;,(.) 

\ cAf) Coif) / \SnJ SNr. 




This form of the NLA amplitude was used in [17] (see also [18]), without account of the 
last two terms in the second line of (43), for the analysis of the total 7*7* cross section. 

Since as we will shortly see the NLA corrections are very large, the choice of the 
representation for the NLA amplitude becomes practically important. In the present 
situation, when an approach to the calculation of the NNLA corrections is not developed 



^The following expression holds actually forn > 1. The coefficient di coincides with rf™^, see Eq. (46). 
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yet, the series representation (38) is, in our opinion, a natural choice. It includes in 
some sense the minimal amount of NNLA contributions; moreover, its form is the closest 
one to the initial goal of the BFKL approach, i.e. to sum selected contributions in the 
perturbative series. 

It is easily seen from Eqs. (38)- (42) that the amplitude is independent in the NLA 
from the choice of energy and strong coupling scales. Indeed, with the required accuracy. 



aM = a.M{l-^^ln(f^^]] (44) 



and therefore terms a" In" ^ s In sq and a" In'^ ^ shinn cancel in (38). 

One can trace the contributions to each d„ coefficient coming from the NLA corrections 
to the BFKL kernel and from the NLA impact factors 

d„ = C + dT" , (45) 
^^ri\n\ — — + -—2 -—In — — + — — / dv f{v)cx{v)c2{i^^ 




\ ^ ' —00 



,n-l 



(46) 



Ci{u) C2(Z/) 

The first coefficient, di, is entirely due to the NLA corrections to the impact factors, 

di^(^r, 4'' = 0. (47) 

Let us note that in the BFKL formalism the NLA contribution to the impact factors 

guarantees not only independence of the amplitude from the energy scale, sq, but it 
contains also a term proportional to In fi^i which is important for the renorm-invariance 
of the predicted results, i.e. the dependence of the amplitude on /ir and Sq is subleading 
to the NLA accuracy. 



3 Numerical results 

In this Section we present some numerical results for the amplitude given in Eq. (38) for 
the Qi = Q2 = Q kinematics, i.e. in the "pure" BFKL regime. The other interesting 
regime, Qi ^ Q2 or vice-versa, where collinear effects could come heavily into the game, 
will not be considered here. We will emphasize in particular the dependence on the 
renormalization scale hr and sq in the NLA result. 

In all the forthcoming figures the quantity on the vertical axis is the L.H.S. of Eq. (38), 
T'ms{A)Q'^ / {s D1D2). In the numerical analysis presented below we truncate the series in 
the R.H.S. of Eq. (38) to n = 20, after having verified that this procedure gives a very 
good approximation of the infinite sum for the Y values Y < 10. We use the two-loop 
running coupling corresponding to the value as{Mz) — 0.12. 
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We have calculated numerically the 6„ and dn coefficients forn/ = 5 and sq = = nj^, 
getting 



17.0664 bi = 34.5920 62 = 40.7609 63 = 33.0618 64 = 20.7467 

65 = 10.5698 ba = 4.54792 67 = 1.69128 bg = 0.554475 

= -3.71087 d2 = -11.3057 ds = -23.3879 d^ = -39.1123 

= -59.207 de = -83.0365 ^7 = -111.151 4 = -143.06 . 



(48) 



In this case contributions to the dn coefficients originating from the NLA corrections to 
the impact factors are 

rf/^P = -3.71087 4°"^ = -8.4361 = -13.1984 df'' = -18-0971 

dfP = -23.0235 CP = -27.9877 4°"^ = -32.9676 4°''' = -37.9618 . ^ ^ 

Thus, comparing (48) and (49), we see that the contribution from the kernel starts to be 
larger than the impact factor one only for n > 4. 

These numbers make visible the effect of the NLA corrections: the dn coefficients are 
negative and increasingly large in absolute values as the perturbative order increases. The 
NLA corrections turn to be very large. In this situation the optimization of perturbative 
expansion, in our case the choice of the renormalization scale fiR and of the energy scale 
So, becomes an important issue. Below we will adopt the principle of minimal sensitivity 
(PMS) [19]. Usually PMS is used to fix the value of the renormalization scale for the 
strong coupling. We suggest to use this principle in a broader sense, requiring in our case 
the minimal sensitivity of the predictions to the change of both the renormalization and 
the energy scales, /ir and sq- 

Since the dependence of results on sq is a feature typical of the BFKL approach and 
is somewhat new for the application of PMS, we will first illustrate the success of PMS 
in this respect on the following QED result known since a long time. In 1937 Racah 
calculated the total cross section for the production of e"'"e~ pairs in the collisions of two 
heavy ions at high energies [20], 



here Zi^2 are the ions charges, mg is electron mass, the ions' four-momenta are pi,2, 

i^ln^i^, (51) 

mim2 

is the energy logarithm and mi^2 are the masses of the ions. The contributions suppressed 
by the power of energy are denoted as C(l/(piP2))- 

The coefficients in front of the subleading logarithms are large and have alternating 
signs 

A = -178/28 = -6.35714 

B = 7^(77r2 + 370) = 15.6817 (52) 
28 
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O/Oo 




Figure 2: c/cio as a function of the energy logarithm / for the cases of exact result of Racah, 
approximated result with Iq = —A/3 (PMS optimal choice) and lo = (kinematical scale 
for energy logarithms). 



C - (348 + yTT^ - 21C(3)^ = -13.8182 . 

To illustrate the application of PMS, imagine that we know only the coefficient A in 
front of the first subleading logarithm. Then using this knowledge we can construct the 
following approximation 

a'^^ = ao ((/ - lo)' + {A + 3/o)(/ - kf) , = ^^^W^' , (53) 

(an analog of NLA in the BFKL approach) where we shift the energy scale introducing 
the parameter Iq. Note that the dependence of the cross section on Zq is subleading in 
the approximation used in Eq. (53). We fix Zq by requiring the minimal sensitivity of 
(53) to the change of this parameter. It is not difficult to find that this procedure gives 
Iq = —A/3 = 2.11905."^ In Fig. 2 we present three curves for cr/cxo as a function of the 
energy logarithm I; the first one was calculated using the exact result of Racah (with all 
subleading logarithms), the other two curves were calculated using (53) with Iq — and 
with the PMS value k = 2.11905. 

We see that the PMS approach gives a very good approximation to the Racah result^. 
On the other hand the procedure with Iq — 0, which means that a kinematical scale for 

^Note that in this example PMS gives the vahie of the parameter for which the correction to the 
lowest approximation, (/ — ^o)^i vanishes. Therefore in this case PMS gives a result which coincides 
with the one given by another alternative approach to optimize the approximation, the fast apparent 
convergence prescription [21]. 

''The negative cross section at Z < 2 is due to the fact that terms subleading in energy, 0(l/(piP2)) 
in (50), are not taken into account. 



11 



2 



4 



6 



8 



Figure 3: 1ms{A)Q'^ / {s D1D2) as a function of Yq at /^r = lOQ. The different curves are 
for Y values of 10, 8, 6, 4 and 3. The photon virtuahty has been fixed to 24 GeV^ 
{rif = 5). 

energy logarithms is used in the approximate formula, makes an awfully bad job for the 
whole I range presented in the figure. 

Returning to our problem, we apply PMS to our case requiring the minimal sensitivity 
of the amplitude (38) to the variation of iir and sq- More precisely, we replace in (38) 
ln(s/so) with Y — Yq, where Y — \n{s/Q'^) and Yq — ln(so/<5^), ^-nd study the dependence 
of the amplitude on 1^. 

The next two figures illustrate the dependence on these parameters for Q'^—24: GeV^ 
and Uf — 5. In Fig. 3 we show the dependence of amplitude on Yq for /ir — lOQ, when 
Y takes the values 10, 8, 6, 4, 3. 

We see that for each Y the amplitude has an extremum in Yq near which it is not 
sensitive to the variation of Iq; or •^o- Our choice of /ir for this figure is motivated by the 
study of ^R dependence. In Fig. 4 we present the /ir dependence for y = 6; the curves 
from above to below are for Yq—3, 2, 1,0. 

Varying /ir and Yq we found for each Y quite large regions in /ir and Iq where the 
amplitude is practically independent on fiR and Yq. We use this value as the NLA result 
for the amplitude at given Y. In Fig. 5 we present the amplitude found in this way as a 
function of Y. The resulting curve is compared with the curve obtained from the LLA 
prediction when the scales are chosen as /ir = lOQ and Yq = 2.2, in order to make the 
LLA curve the closest possible (of course it is not an exact statement) to the NLA one in 
the given interval of Y. The two horizontal lines in Fig. 5 are the Born (2-gluon exchange) 
predictions calculated for /ir — Q and /ir — lOQ. 
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Figure 4: Ims{A)Q'^ / {s D1D2) as a function of at Y=Q. The different curves are, from 
above to below, for Yq values of 3, 2, 1 and 0. The photon virtuality has been fixed to 



Similar procedure was applied to a lower value of the photon virtuality, Q^=b GeV^ 
and Uf — 4, where we again found that the choice of scales jiR — lOQ and Yq — 2.2 makes 
the LLA amphtude almost the closest to the NLA curve. The results are presented in 



Note that in calculating the NLA amphtude for Figs. 5, 6 we use optimal values of jiR 
and Yq found to be a bit different for different Y values from jiji — lOQ and Yq — 2.2. 

We stress that one should take with care BFKL predictions for small values of Y, since 
in this region the contributions suppressed by powers of the energy should be taken into 
account. At the lowest order in ag such contributions are given by diagrams with quark 
exchange in the t-channel and are proportional in our case to aEM(^sfv/Q'^- At higher 
orders power suppressed contributions contain double logarithms, terms ~ In^"^ s, which 
can lead to a significant enhancement. Such contributions were recently studied for the 
total cross section of 7*7* interactions [22] . 

If the NLA (and LLA) curves in Figs. 5,6 are compared with the Born (2-gluon ex- 
change) results, one can conclude that the summation of BFKL series gives negative 
contribution to the Born result for y < 6 if one chooses for the scale of the strong cou- 
pling in the Born amplitude the value given by the kinematics, iir — Q. We believe 
that our calculations show that one should at least accept with some caution the results 
obtained in the Born approximation, since they do not give necessarily an estimate of the 
observable from below. 

Another important lesson from our calculation is the very large scale for ag (and 
therefore the small itself) we obtain using PMS. It appears to be much bigger than the 
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Figure 5: Ims{A)Q'^ / (s D1D2) as a function of Y for optimal choice of the energy param- 
eters Yq and (curve labeled by "NLA"). The other curves represent the LLA result for 
Yq = 2.2 and iir = lOQ and the Born (2-gluon exchange) hmit for Hr = Q and hr = lOQ. 
The photon virtuality has been fixed to 24 GeV^ (n/ = 5). 



kinematical scale and looks unnatural since there is no other scale for transverse momenta 
in the problem at question except Q. Moreover one can guess that at higher orders the 
typical transverse momenta are even smaller than Q since they " are shared" in the many- 
loop integrals and the strong coupling grows in the infrared. In our opinion the large 
values of /ir we found is not an indication of the appearance of a new scale, but is rather 
a manifestation of the nature of the BFKL series. The fact is that NLA corrections are 
large and then, necessarily, since the exact amplitude should be renorm- and energy scale 
invariant, the NNLA terms should be large and of the opposite sign with respect to the 
NLA. We guess that if the NNLA corrections were known and we would apply PMS to 
the amplitude constructed as LLA -|- NLA-corrections + NNLA-corrections, we would 
obtain in such calculation more natural values of /ir. 

We conclude this Section with a comment on the possible implications of our results 
for mesons electroproduction to the phenomenologically more important case of the 7*7* 
total cross section. By numerical inspection we have found that the ratios bn/bg we got for 
the meson case agree forn = 1-rlO at 1-^2% accuracy level with the analogous ratios for 
the longitudinal photon case and at 3.54-30% accuracy level with those for the transverse 
photon case. Should this similar behavior persist also in the NLA, our predictions could 
be easily translated to estimates of the 7*7* total cross section. 
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Figure 6: The same as in Fig. 5 for photon virtuaUty Q"^ fixed to 5 GeV^ (n/ = 4). 



We have determined the amphtude for the forward transition from two virtual photons to 
two light vector mesons in the Regge limit of QCD with next-to-leading order accuracy. 
This amplitude is the first one ever written in the next-to-leading approximation for a 
collision process between strongly interacting colorless particles. It is given as an integral 
over the v parameter, which labels the eigenvalues of the leading order forward BFKL 
kernel in the singlet color representation. This form is suitable for numerical evaluations. 
The result obtained is independent on the energy scale Sq, and on the renormalization 
scale [iR within the next-to-leading approximation. 

Using a series representation of the amplitude which includes the dependence on the 
energy scale and on the renormalization scale at subleading level, we performed a numeri- 
cal analysis in the kinematics when the two colliding photons have the same virtuality, i.e. 
in the "pure" BFKL regime. We have found that the next-to-leading order corrections 
coming from the kernel and from the virtual photon to light vector meson impact factors 
are both large and of opposite sign with respect to the leading order contribution. 

An optimization procedure, based on the principle of minimal sensitivity method, has 
proved to work nicely and has lead to stable results in the considered energy interval, 
which allows us to predict the energy behavior of the forward amplitude. The procedure 
consists in evaluating the amplitude at values of the energy parameters for which it is the 
least sensitive to variations of them. We have found that there are wide regions of values 
of So and \xr where the amplitude remains almost fiat. 

The optimal choices of sq and \ir are much larger that the kinematical scales of the 
problem. More than being the indication of appearance of another scale in the problem. 



4 Conclusions 
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this could be related to the nature of the BFKL series. The renorm- and energy scale 
invariance, together with the large next-to-leading approximation corrections, call for 
large next-to-next-to-leading order corrections, which are most probably mimicked by 
unnatural optimal values for sq and /ir. 
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